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1. Introduction 

The European Twisted Mass (ETM) collaboration is pursuing its effort to generate Nf = 2 + 
1 + 1 configurations of twisted mass fermions. This discretization of QCD has the advantage of 
providing automatic 0{a) improvement of all the observables once tuned at maximal twist [|l|, |2|]. 
In this contribution we are interested in computing the masses of the baryon octet and decuplet. 
The inclusion of the strange and charm quarks in the sea offers the possibiUty to investigate several 
effects in the baryon sector. 

First, given the large value of the charm quark masss, potentially large cutoff effects could 
appear in observables believed to be insensitive to the charm degrees of freedom. We perform a 
preliminary analysis of the discretization effects using three lattice spacings. The non degener- 
ate twisted mass fermion doublet used in the heavy sector do not conserve the strange and charm 
quantum numbers at finite lattice spacing. The effects of this symmetry breaking have already been 
investigated in [^] to determine the kaon and D meson masses. In order to study the spectrum of 
strange baryons, two approaches can nevertheless be followed: the mixed action setup using an ac- 
tion that conserves strangeness in the valence or the unitary setup using the sea action. Both meth- 
ods have been followed to investigate unitarity violations. Finally, by comparing present results 
with our previous work using Nf = 2 configurations in the partially quenched approximation [Q], 
we investigate the effect of the quenching of the strange quark. 

2. Lattice setup 

In the gauge sector we use the Iwasaki gauge action |^. Using the same notations as in 
Ref. [^, the light doublet action reads : 

Si = a^Y.{zi{x)[D[U]+mo + in,r5r3]Ziix)} , (2.1) 

X 

while the action in the heavy non degenerate doublet is : 

Sh = «'^X{^''W[^[^]+'«0 + iAtCT75Tl+jU5'!^3]Z/!W} , (2-2) 

X 

The input parameters of the calculation (L/a, jS,/^/,/!^,/^^) are collected in Table [I[ They 
span a pion mass range from approximately 270 to 500 MeV and mpsL > 3.4 ( the lowest bound 
holds for ensembles B25.32 and A40.24). Estimation of the lattice spacing from Ref. ||^ gives 
fl'j3=i.90 = 0.0861(4) fm and fl:jg=i.95 = 0.0778(4) fm. A preleminary estimate of the scale at 
j3 = 2.10 gives ap^2.w ~ 0-060 fm. 

The masses are extracted from two-points correlators using the same interpolating fields as in 
our previous Nf = 2 study Errors are estimated using the bootstrap method. 

To improve the overlap of the interpolating fields with the ground states we apply Gaussian 
smearing to each quark fields defining q{x,t) = Y.y ^ {x,y;U {t))q{y ,t) using the gauge invariant 
smearing function 

F{x,y;U{t)) = {\ + aH)"{x,r,U{t)), (2.3) 
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Ensemble j3 a/i/ a/ifj aiJ.§ {L/aYxT/a 

A30.32 L90 0.0030 0.150 0.190 32^ x 64 

A40.32 0.0040 32^ x 64 

A40.24 0.0040 24^ x 48 

A50.32 0.0050 32^ x 64 

A60.24 0.0060 24^ x 48 

A80.24 0.0080 24^ x 48 

A100.24 0.0100 24^ x 48 

B25.32 h95 0.0025 0.135 0.170 32^ x 64 

B35.32 0.0035 32^ x 64 

B55.32 0.0055 32^ x 64 

B75.32 0.0075 32^ x 64 

B85.24 0.0085 24^ x 48 

D20.48 110 0.0020 0.120 0.1385 48^ x 96 

D30.48 0.0030 48^ x 96 



Table 1: Summary of the A^f = 2 + 1 + 1 ensembles used in this work. 

constructed from the hopping matrix, H{x,y;U{t)) = Y]=[(^i{'^,l)Sx,y-i + Uj {x — i,t)dx,y+i^ ■ 

Furthermore we apply APE smearing to the spatial links that enter the hopping matrix. The param- 
eters for the Gaussian and APE smearing are given in Table ^. 



p 


^ Gaussian 


O^Gaussian 


A^APE 


ccape 


1.90 


50 


4.0 


20 


0.5 


1.95 


50 


4.0 


20 


0.5 


2.10 


70 


4.0 


20 


0.5 



Table 2: Smearing parameters used in this work. 

Two approaches are followed to investigate strange baryons using our Nf = 2 + 1 + 1 simula- 
tions. 

The first one consists in using a mixed action in order to avoid the mixing of the strange and 
charm quantum numbers. As in our Nf = 2 study, we choose Osterwalder-Seiler fermions in the 
valence to simulate the strange quark. They are defined by the following action: 

= £ {D[U]+ mo + i^sYs ) six) (2.4) 

X 

The parameter aj^s has been tuned using the kaon mass calculated in the unitary theory on each 
gauge ensembles by imposing : 

{ajis,^ ,a}Xua}Xa ,alXs) = ni'^^^'^ {^,apLi,alXa,a}Xs) (2.5) 

In practice one performs the measurement of the baryon masses for several values of the bare 
strange quark mass apis and interpolate linearly to the matched quark mass. All the results pre- 
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sented in the following are interpolated at the matched strange quark mass. 

The second method is using the unitary theory as defined in Eq. (p^. The measurements of 
the kaon and D meson masses in the unitary setup have been extensively discussed in Ref. [^]. In 
particular, one can show that performing a rotation to the physical basis assuming all the renormal- 
ization are equal to one (the so called pseudo physical basis) leads to a propagator that simulates 
a strange quark for large Euclidean time. We use this strategy to extract the octet and decuplet of 
strange baryons. 

Note that the two approaches differ by ^{a^) effects and that the difference in the results 
between the two approaches is an unitarity violation. 



3. Results in the light sector 

We show in Fig. |l|, our results for the nucleon mass in ro units as a function of (romps )^. The 
physical point is shown using ro = 0.447(5) fm from Ref. [^. In this plot as well as in all the 
following ones, a systematic error on the plateau determination has been added quadratically to the 
statistical error. The results obtained are consistent, indicating that finite volume and lattice spacing 
effects are small. 

Another important issue raised by twisted mass fermions, is the isospin symmetry breaking at 
finite lattice spacing. In the Nf = 2 case we have shown in Ref. [0] that in the A sector, isospin 
symmetry breaking manifests itself as a mass spliting between A++'^ and A+'° states and was 
found to be compatible with zero. The same conclusion holds for A/y^ = 2 + 1 + 1 results as shown 
in Fig. ^ where we display the relative mass difference between the two states for three j8 values 
as a function of romps. 
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Figure 1: Nucleon mass in units of ro as a function 
of (roOTps)^. Physical point is indicated by a black 
star and the scale has been set using ro(j3 ~ 1.95) = 
0.447(5) fm from [|l 



Figure 2: Difference of the A++' and a A+' masses 
divide by their average as a function of romps 



Given that we have results for three lattice spacing, it is worthwhile to perform a combined fit 



4 



Low lying baryon spectrum with Nf = 2+1 + 1 dynamical twisted quarks 



Vincent Drach 



of our results. To this aim we used the fitting function 



rom^v = rom^N^ + ^{rmpsf - jMni^omps? +C{-)^ (3.1) 



where the parameter is the nucleon mass in the chiral limit, cj^ is related to the nucleon sigma 
term, gA is nucleon axial coupling, fj^ the pion decay constant and C is a coefficient that controls 
the discretization effects. In our fits, the values of gA and fn are fixed to their experimental values 
i.e gA = 1-26 andfjc = 0.1307 GeV. This parametrization comes from an &{p^) expansion in heavy 
baryon chiral perturbation theory (see Ref. [Q] and reference therein for details). In order to take 
into account possible discretization effects a term proportional to has been added in Eq. (BTl|). 
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Figure 3: Nucleon mass in ro units as a function of 
(romps )'^. The red curve comes from the combined fit 
using Eq. (3. 1) of all the Nf = 2 + \ \ data, setting 
fl = 0. The blue dotted lines correspond to the same 
results with Nf=2. The physical point is not included 
in the fit. 



Figure 4: Nucleon mass as a function of (a/ro)^ for a 
fixed pion mass reference. The green dots are the nu- 
cleon masses interpolated at romps = 0.73. The blue 
dotted line is obtained from our combined fit. The 
continuum exttapolation is shown by a black triangle. 



We show in Fig. ^ the nucleon mass as a function of (romps)^. The lattice results were 
fitted using Eq. (3.1) and the best parametrization extrapolated to the continuum (setting a = ) is 



denoted by a red curve. The fit has a ^ = 1 1 .2/10. The finite lattice spacing data are compatible 
with our estimation of the continuum behaviour, showing that the lattice spacing dependence is 
small. In order to assess for possible quenching effects of the strange quark we show with a blue 
dotted line the same fit performed using our Nf = 2 data at three lattice spacings (xlo f = 5.5 /7 ). 
The two fits agree with each other, indicating that, within errors, the nucleon mass is not sensitive 
to the presence of the strange quark in the sea. 

The scaling behaviour of the nucleon mass, for a fixed pion mass reference, is shown in Fig. ^. 
The value romps = 0.73 is in a region where we have data for the three lattice spacings. The 
finite lattice spacing points have been interpolated from the data to match the reference pion mass. 
The blue dotted line indicates the best fit parametrization of our results. Our estimation of the 
coefficient C is compatible with zero, confirming that discretization effects are smaller than our 
statistical errors. 
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4. Results in the strange sector 



In our previous study on the strange baryons [Q] we have shown that the isospin breakingin 
was manifests itself in the S sector. We show in Fig. ^ the mass spUtting for the S as a function of 
romps for the three lattice spacings. Most of our results are now compatible with zero. 
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Figure 5: Mass splitting of the E as a function of romps for the three lattice spacings. 



In Fig. 1^ we show the A mass, computed in the mixed action approach explained in in units 
of ro as a function of (rowips)^ for the two finest lattice spacings. We also show our Nf = 2 results 
obtained at j8 = 3.9. The results agree within the errors. A quantitative comparison is made difficult 
by the fact that the sea strange quark mass at j8 = 1.95 and j8 = 2.10 are not matched. Note also 
that the valence strange quark mass used in our Nf = 2 computation has not been estimated. With 
the existing data, we are unable to obtain quantitative estimates of the discretisation effects in our 
Nf = 2 + 1 + 1 results or to highlight some "quenching effect" with respect to the Nf = 2 case. 
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Figure 6: A mass as a function of (romps )^ in the 
mixed action approach. For comparison Nf = 2 data 
are shown with black dot. The physical point is in- 
dicated by a black star. The scale was fixed using ro 
from Ref . [@] 



Figure 7: S mass as a function of (romps)^ both in 
the unitary and mixed action approach. For compari- 
son Nf = 2 data are shown with black dot. The phys- 
ical point is indicated by a black star. We used the 
value of ro from Ref. [S] to set the scale. 
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In Fig. ^ besides the mass of the S as a function of (romps)^ obtained in the mixed action 
approach and in the Nf = 2 case, we show the results obtained in the unitary setup. We observe 
that there is a good agreement except for ensemble B55.32 where the unitary theory gives a lower 
value. The behaviour at this particular value is under study. 

5. Conclusion 

We have shown preleminary results concerning the ground state baryon spectroscopy using 
Nf = 2 + I + I simulations. The nucleon mass shows small discretisation effects and a chiral 
behaviour similar to the one obtained in our Nf = 2 study. To our present accuracy, the isospin 
breaking is compatible with zero both in the A and E case. We briefly described the two strategies 
followed to determine masses in the strange sector, namely the so called "mixed" and "unitary" 
cases. The results agree within errors, but it is stressed that with the present amount of informations 
quantitative statements on discretisation or quenching effects would be premature. In all case the 
results obtained with the mixed approach are in agreement with those in the unitary case, with one 
exception in the case of the S where further investigation is under way. 
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